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Abstract 

In this work we investigate structural properties of native states of a 
simple model for short flexible homopolymers, where the steric influence 
of monomeric side chains is effectively introduced by a thickness con- 
straint. This geometric constraint is implemented through the concept 
of the global radius of curvature and affects the conformational topology 
of ground-state structures. A systematic analysis allows for a thickness- 
dependent classification of the dominant ground-state topologies. It turns 
out that helical structures, strands, rings, and coils are natural, intrinsic 
geometries of such tubelike objects. 

PACS: 05.10.-a: Computational methods in statistical physics and non- 
linear dynamics; 87.15.A-: Theory, modeling, and computer simulation; 
87.15.Cc: Folding: thermodynamics, statistical mechanics, models, and 
pathways 

1 Introduction 

The structural properties of macromolecules composed of covalently bound 
atomic complexes are of major interest as the functionality of these molecules, 
e.g., biopolymers, strongly depends on the formation of stable ground-state 
conformations with several substructures. One can resolve these structures by 
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means of several experimental techniques like molecular microscopy, NMR or 
X-ray analyses. 

As these techniques are relatively costly and can often hardly be general- 
ized, the structural behavior of polymers and its modeling got into the focus 
of computer simulations, too. In early approaches, polymers were modeled 
as topologically one-dimensional strongly coarse-grained bead-spring systems, 
which were treated by means of Monte Carlo and molecular dynamics com- 
puter simulations JXJ [21 E] - Later, all-atom peptide simulations were performed 
with the intention to study real biopolymers, in particular their native states 
and the folding pathway to native states [H [5]. As this is still a great chal- 
lenge and moreover restricted to comparatively short objects, the interest in the 
earlier simple models continues until today. There is still a great variety of one- 
dimensional stringlike models of flexible polymers under current investigation 
(see, e.g. H E]). 

Biopolymers do have side chains, however, which are responsible for the vari- 
ety of forms and functions. As these side chains involve strict steric constraints 
for the structure of the polymer, the central question is therefore: What de- 
gree of abstraction (coarse-graining) is reasonable to treat certain features of 
real (bio)polymers or, the other way around, which features can we reliably 
study with a certain degree of abstraction? One can, for example, understand 
certain universal properties of the well-known coil-globule transition by study- 
ing the simplest stringlike models [TUl [H] , whereas the formation of secondary 
structures will generally not be answered satisfactorily^ 

We are going to approach this problem in the present work by studying in de- 
tail some kind of "intermediate" model. It is derived from a linelike model where 
additionally the steric influence of monomeric side chains in real biopolymers 
will be effectively introduced, without taking into account further microscopical 
details, by a geometric thickness constraint, extending the model to a three- 
dimensional tube [TBI [TTl [18] . An elegant possibility to implement the thickness 
constraint is provided by the concept of the global radius of curvature which 
"provides a concise characterization of the thickness of a curve" [IS]. We will 
see, that we already get quite interesting qualitative and quantitative answers 
to the question of secondary structure formation under these conditions]! 

The structure of this paper is as follows: In Sect. 12.11 we first introduce the 
model we use and shortly describe our methods, and in Sect. 12.21 we explain in 
more detail, how the thickness is implemented in our model. After introducing 
briefly some observablcs in Sect. 13. 1[ we study in Sect. 13.21 the ground states 
of the model in great detail, depending on the thickness with inter-monomer 
interaction parameters kept fixed. Furthermore, in Sect. 13.31 we look at special 



We refer here mainly to ground states of models for flexible polymers with a single length 
scale. Of course, there are studies of simple polymer models dealing with the problem of 
secondary structure formation, where helices were found for example as (long-living) transient 
states during the folding process, as ground states for stiff polymers, or at special ratios of 
different length scales in the system 1121 1131 1141 115| . 

2 There are recent, interesting studies of a tube model, showing that helices can form en- 
tropically driven. In these studies, however, solvent particles and solvent effects are explicitly 
introduced and play an important role 1201 1211 1221 . 
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problems like crystallization into regular lattice structures and the appearance 
of biologically relevant structures depending on a variable interaction length 
scale. The paper concludes in Sect. 0] with a summary of our main findings. 



2 Technical details 
2.1 Model and Methods 

To model flexible homopolymers we use a thick, tubclikc off-lattice chain with 
fixed bond length and pure Lennard- Jones (LJ) interaction between all, except 
for the neighboring, monomers: 



where = |xj — Xj| is the distance between two monomers at positions Xj 
and Xj. Note that VLj^ij) = for = a and that the potential minimum 
Vhj{rij) = -1 is located at r^ in = 2 1 / 6 cr. The bond length r i)i+ i is set to 1. 

We investigate the model by means of Monte Carlo computer simulations 
and numerical methods. For the ground-state search, one may use for exam- 
ple the multicanonical method [231 EH E5] , the efficient random walk algorithm 
introduced by Wang and Landau [26] [27] or parallel tempering [28l [29] . For 
this purpose, all methods work nearly equally well as one does not have to 
care about the quality of the sampling of the whole configuration space and the 
performance depends mainly on the parameters of the methods. Additionally, 
we use standard conjugate gradient methods to refine the results [30]. For the 
methodologically more challenging task of studying the thermodynamic behav- 
ior [2I1[32], w e use the multicanonical method, calculating the multicanonical 
weights with the help of the Wang-Landau algorithm [25] . Chains with different 
thicknesses have always been simulated independently from each other, not at 
least to avoid uncontrollable correlations. 

2.2 Thickness and Global Radius of Curvature 

What precisely motivates us to simulate coarse-grained homopolymer models 
with geometric constraints like "thickness"? And, how is this realized in the 
model? Polymers in biology are not thin strings. Amino acids, and thus pro- 
teins, do have rather extended side chains sterically avoiding each other. It 
might therefore be useful to introduce a constraint that mimics this volume 
exclusion. Furthermore, it has been shown [33] that tube models for polymers 
allow for the formation of stable regions of biologically relevant (sub) structures 
like helices or sheets, in contrast to simpler linelike polymer models [M]. In 
addition, the introduction of such a geometric constraint restricting the confor- 
mational space might even lead to some technical advantages for finding minimal 
energy conformations in sophisticated stochastic ground-state searches [32] . 





(1) 
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To define the self-avoiding tube, we use the concept of the global radius of 
curvature r gc , which is a concise characterization of the thickness of a curve [191 
IM] . The global radius of curvature is defined as the minimal radius of all 
circumcircles r c of any three monomers in the chain: 

r gc (X) := min{r c (x. t , x jt x fe ), V 1 < i < j < k < N} . (2) 

The circumradius r c , i.e., the radius of curvature, of three points located at 
positions x^, Xj and x& can be calculated as 

rcCxi.x^Xfc) = — -j J - , (3 

where (x, , x 3 ■, x^ ) is the area of the triangle spanned by the three points at 
Xj, Xj and x&. Note, that the model uses a three-point interaction between 
monomers, in contrast to two-point interactions typically used in hard-sphere 
models to incorporate volume exclusion effects. Given a polymer conformation 
X = (xi, . . . , x;v) with N monomers, we then define as its thickness d twice the 
global radius of curvature r gc (X): <f(X) = 2r gc (X). 

The other way around, given a thickness constraint p such that d > 2p, one 
can construct an excluded volume depending on p around two monomers, which 
is "forbidden" for any other monomer. A polymer conformation then complies 
with the thickness constraint if any other monomer resides outside these circles. 
The partition function reads in this case: 

Z p = jvxe(r ec (X)-p)c- E ^/ T , (4) 

where Q(z) is the Heaviside function and 

E(X)= J2 V u {r t3 ) (5) 

is the energy of the conformation X. 

Intuitive illustrations of this approach can be found in |18j . where it has 
been chosen for the analysis of ringlike tube polymers. It is used as well in [3T] 
[32] where the pseudophases of secondary structures for tubelike polymers are 
investigated. As a remark: Even though some universal properties (such as, 
e.g., critical exponents) do not depend on the exact definition of the thickness, 
the model may behave differently by imposing the hard constraint d = 2p, i.e., 
by replacing 0(r gc (X) — p) with 5(r gc (X) — p) in the partition function (jU [TS] . 



3 Results and Discussion 

In the following, we first analyze in detail the ground states and ground-state 
regions of tubelike polymers for chain lengths 8 < N < 13 with Lcnnard- Jones 
parameter a = 1, where Vlj = for covalcntly bound monomers. Furthermore, 
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we investigate the influence of the tube thickness in connection with the length 
scale a of the nonbonded LJ interaction on the formation of different classes 
of structures. In this analysis, particular emphasis will be dedicated to sec- 
ondary structures such as helices and strands, being most relevant for segments 
of biomolccules. 

3.1 Conformational Classification Observables 

We characterize and identify conformations by their energy E(X), but naturally 
also by geometrical properties such as the end-to-end distance r cn d or radius of 
gyration r gyr . However, as this turns out to be not sufficient to distinguish 
between all conformations, we also take into account local radii of curvature 
r \c,i : " r c( x i> x i+ij x i+2); which are related to the bending angles between two 
bonds •& via d = l/?*i c + 0(rj~ 2 ) for small as well as torsion angles <fi G (— 7r, n]. 
We speak of a Ko-conformation, if the chain has a constant local curvature 
at all monomer positions, i.e., ri Cl , = const , Vi. Analogously, we call a structure 
with constant torsion angles a ro-conformation |35j . For example, a prominent 
structure with both, kq- and ro-property, is the perfect a-helix. A conformation 
is called "closed" if the distance between both ends of the chain, i.e., r on( j, resides 
in the close vicinity of the Lennard- Jones minimum, whereas a "symmetric" 
structure exhibits a symmetry of the torsion angles with respect to the center 
of the chain. Nice examples for "closed" kq conformations are twisted circles 
of constant curvature (in German, so-called "windschicfe Kreise" [36]). Finally, 
in "flat" conformations, the backbone has an almost two-dimensional, planar 
structure, where all torsion angles converge to 0, i.e., <f& 2 ^' wncrc 

:= min(|0|,7T- |0|) . (6) 

3.2 Thickness Dependent Ground-State Properties 
of Tubelike Polymers 

For our comprising study of the ground states of the model, we first set the 
Lennard- Jones parameter to a = 1. In this case, the "natural thickness" r™ n , 
i.e., the global radius of curvature of the ground state of a flexible LJ poly- 
mer without thickness constraint, is about half the interaction length r™ ln /2 = 
2-5/6 0.56, which thus sets a reasonable bound for the thickness constraint^ 
Below this value, the thickness constraint docs not influence ground-state prop- 
erties at all. Thus, in the following, we only consider tube polymers with 
P > rp n - 

Figure Q] shows the energies of the ground states for various chain lengths 
in dependence of p in comparison to the energy of the corresponding calcu- 
lated space-filling (perfect) a-helix. This helix is defined as the helix (x,y,z) = 
(r sin </>, r cos <fi, p4>/2iv) with a pitch p such that the surface of the tube has a self- 
contact at the cylinder with radius r and the radius r is minimal under the thick- 
ness constraint. In other words, the optimally packed, i.e., space-filling, a-helix 

3 Actually, due to the discrete nature of the bead-stick polymer model, we measure a 
"natural thickness" > 0.59 for all polymers considered. 
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Figure 1: Energies of ground states depending on the thickness constraint p. 
Dashed lines show for comparison the energy for exact a-helices. 



corresponds to the transition between the two qualitatively different regimes of 
p/r > c* and p/r<l |37j . The computation of this helix is non-trivial, as the 
critical ratio c* depends on the discretization level and ranges from c* « 2 for 
p « 0.7 to c* = 2.512 for the continuous case, which is equivalent to p — > oo. An 
interesting and detailed discussion of compact helix formation and the critical 
ratio c* can be found also in [2"T] . 

For two exemplified chain lengths, N = 8 and N = 10, we have plotted in 
Fig. [2] the derivatives dE/dp in order to emphasize regions of structural activ- 
ity. In these regions, where the derivative exhibits peaks, noticeable qualitative 
conformational transitions occur. To describe and understand these different 




Figure 2: Numerical derivatives of the energies in Fig. [T] with respect to p for 
N = 8 and 10. 
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Figure 3: Ground-state conformations for selected thickness parameters p for 
N = 8, 9, 10 and 13 (from top to bottom). The second row for TV = 8 shows the 
same ground states with appropriate thickness, to give a better idea of the "real" 
objects we are investigating. The second rows for TV = 9 and N = 13 show an 
alternative view of the same configuration. For reasons of better visibility, the 
thickness is not shown in the proper scale (except for the second row). 

classes of ground-state conformations, we visualize in Fig. [3] significant struc- 
tures and plot in Fig. @] the contact maps, where a contact is counted, if the 
distance between two monomers r,j < r™ ln + e. We set e « 0.2 but, of course, 
the contact maps do not depend on minor variations of eQ In addition, in 
Figs. [5] and [5] end-to-end distances and mean torsional angles discussed above 
are shown. Based on this data, we can classify the generic behavior by intro- 
ducing three general regions: thin, intermediate, and thick tubes. Let us now 
look at these regions in more detail. 



4 Furthermore, due to the small size of the systems, the contact maps do not become more 
meaningful by scaling e with the thickness in some way, instead of keeping e constant. 
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Figure 4: Contact maps of the conformations shown in Fig. [3] The axes show the 
monomer numbers i and j, the entries of the matrix are set (•), if two monomers 
i and j are in contact with each other, i.e., if the distance between them in the 
three-dimensional structure is shorter than a certain threshold value. This value 
is here slightly larger than the minimum distance of the LJ potential rf^ ln . 



3.2.1 Thin tubes 

The thin-tubes region is, besides some singular points, dominated by helical 
and helicallike conformations. We call a conformation "helical", if (j> = 0^/2, 
where x is the average along the chain, x = (1 /iV) Xi, i.e., if all torsion 
angles lie in a range . . . ±ir/2 (cp. Eq. © and Fig. [5] we do not distinguish 
between right- and left-handed helices, but the sign must not change within the 
conformation). Furthermore, the entries in the contact map lie precisely parallel 
to the diagonal of the matrix in these clear indication for helix structures 

(see, for example, N = 8; p = 0.685 or N = 9; p = 0.73 in Fig. gj}. "Helicallike" 
conformations share some properties with helical structures, e.g., they exhibit 
a large, slightly increasing end-to-end distance (cp. Fig. [5]) with increasing 
thickness, but the torsion-angle criterion above may be violated (typically in a 
periodical manner) and the contact-map entries do not form an exact parallel, 
but a line roughly parallel to the diagonal (for example, at N = 8; p = 0.6 or 
N = 9;p = 0.78). Generally, we find three interesting effects looking at the 
contact maps in Fig. HI which have been mentioned above or will be discussed 
later again: First we see, that polymer chains without thickness constraint 
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Figure 5: End-to-end distances of ground-state conformations depending on the 
thickness constraint p. The dashed line for N = 9 shows for comparison the 
end-to-end distance for the exact a-helix. 
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Figure 6: Mean torsional angles of ground-state conformations depending on 
the thickness constraint p. Shown are absolute mean values of <p and <\>kIi := 
min(|0|,7r — |0|). The short-dashed line for N = 9 shows the behavior for the 
exact a-helix. 



(see the maps for TV = 8 and TV = 9 with p = 0.6) do not have a pronounced 
structure. Just by increasing the thickness a bit, clear helical structures emerge, 
indicated by straight "lines" parallel to the diagonal of the map (cp. TV = 8 
and TV = 9 with p = 0.685(5)). Secondly, by increasing the thickness further, 
we see, for example for TV = 9, that this parallel "lines" moves away from the 
diagonal, i.e., the helical conformations are "untwisting". Finally, looking at the 
contact maps for TV — 13, we see that tertiary effects come into play, indicated 
by "disrupting" vertical "lines" , which is typically an indication for sheetlike 
structures. 

Remarkably within certain intervals (TV = 8: 0.63 < p < 0.688; TV = 9: 
0.673 < p < 0.6855), the ground-state conformations expand with increasing 
thickness to a perfect space-filling helix with kq- and ro-property, i.e., an a- 
helix with constant bond- and torsion angles {(jj 2 — 4>(j> an d t 2 \c — He He vanish^) . 
The comparison of measured observables with the data for the exact a-helix is 
emphasized in the insets of Fig. Q] and furthermore exemplarily shown in Figs. [5] 
and [51 We will resume the discussion on this fact in Sect. 13.3.21 below. 

As a remark: What is the motivation to call these conformations a-helical 
in imitation of the real biological a-hclix? In natural proteins, a-helices possess 
about 3.6 amino acids per helix turn |38| and have mainly constant bond and 
torsion angles. If we construct a perfect space filling helix with exactly 3.6 
monomers per turn, we find that it has a global radius of curvature of r gc w 0.69. 
Or, the other way around in the region 0.6845 < p < 0.688 (example for TV = 8), 
one counts 3.576 . . . 3.596 monomers per turn, which is in very good agreement 
with natural a-hcliccs. We thus see the first biological relevant structure realized 
by the simplest model with just Lennard-Jones interaction and thickness but 
without any conformational assumptions or additional input. 

An above mentioned singular point is located in the vicinity of the perfect 
helices at p m l/y/2 « 0.71, where ground states attempt to crystallize in a 
regular simple cubic (sc) lattice structure. We find for example for TV = 8 at 
p = 0.73 a Ko-conformation almost fitting the sc lattice (elsewhere called "simple 
cubic lattice helix" [33]), which then untwists with increasing thickness. We see 
the same tendency for longer chains as well (see Fig. [3] for visualizations and 
Scct. l3.3.2"l for further discussion). Note that a perfect cube will not be a ground 
state at any thickness, as the Lennard- Jones interaction length scale is larger 
than the bond length. If we reconfigure the potential such that its minimum 
value equals the bond length, i.e. set r™ ln = 1, we find indeed that the ground- 
state conformations fit exactly into the simple cubic lattice (i.e. are exact cubes 
for adequate monomer numbers) up to lengths of ~ 30 monomers. We will show 
this in more detail in Sect. 13.3.11 

At larger thickness we observe in Fig. [3] extended helicallike conformations, 
which may overlap due to the shortness of the chains only at the end bonds. 

5 A remark on the precision of the simulation: The values of <j> 2 — 4><j> and r 2 i c — f\ c f\ c 
become even with the stochastic methods smaller than 10 -8 at this point, i.e., the difference 
between any two torsion angles, for example, in the chain is already less than 1.5 X 10 — 4 7r. 
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3.2.2 Intermediate tubes 



In the interval 0.9 < p < 1.0, we observe an abrupt switch to almost flat 
(cp. Fig. and mostly closed (cp. Fig. [5]) conformations. One finds bended 
double-rings, hairpins, and even conformations that are "crystallized" on a two- 
dimensional honeycomb lattice (cp. Fig. [31 N = 8, 13, p w 1.02). These curves 
are, of course, Ko-curves as well and have apparent similarities to /3-sheets known 
from secondary structures of biopolymers. We find in some small regions com- 
petitions between mesomeric structures, i.e. structures with the same monomer 
positions but different bond distributions (see, for example, Fig. [3J N = 10, 
pwl.l). 

3.2.3 Thick tubes 

At p ss 1.1, the ground-state conformation is (again) "closed" for all chain 
lengths. Here begins the region of the twisted circles of constant curvature 
( "windschiefe Kreise") [3(3 [35]. With increasing thickness, the rings become 
more and more flat until they reach the two-dimensional ring at p w N/2tt, 
which is again a kq- and ro-curve. Increasing thickness just pushes apart the 
ring, what can clearly be seen in the end-to-end distance and the torsion angles 
(see Figs. [5] and [6]). For the somehow pathological case of p — > oo one would 
reach the limit of stiff rods. 

For purposes of illustration we display two examples from the class "wind- 
schiefe Kreise" . The first kind consists out of 4 half circles, which form a closed 
three-dimensional curve. The left side of Fig. [7] displays a N = 32 chain, which 
is a ground state of the theory and has been obtained from simulations at 
r gc = 2.562915 and r™ ln = 1.6. As can be seen each of the half-circles consists 
of eight monomers. The second kind of "windschiefe Kreise" consists out of four 
helix sections, that arc joined together in such a way that the resulting curve is 
closed again. The right side of Fig. [7] displays such a N = 32 chain, which was 
obtained from simulations at r KC = 3.624510 and r'" m = 1.6. 




Figure 7: Two examples for twisted circles of constant curvature r c m 2.6 (left) 
and r c w 3.6 (right) with = 32 monomers. See text for details. 
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3.2.4 General remarks 

It is not surprising, that the situation becomes more complex with increasing 
chain length. At least some of the described "nice-looking phases" above are 
artificial in the sense, that they occur at exactly one short length, or are favored 
just by that very short length, respectively. We see, for example, for N = 10 
and N = 13 no exact (a-)helices anymore, it rather seems that at these lengths 
"tertiary" effects already play a role in the sense, that two small secondary 
structures are formed which are then arranged "side by side" . An indication for 
this trend may be that conformations with low thickness are often "symmetric"!!, 
i.e., the conformations get buckled and turn back at some point (generally in the 
middle). See for example the helical region for N = 13 in Fig. [3] Anyhow, the 
helical structures being present for shorter chains indeed exist "very close" to 
the ground states, i.e., with a slightly higher energy. Two of these conformations 
are depicted in Fig. [5] 

We will get further convincing arguments for this classification scheme by 
investigating the thermodynamic behavior of these polymers in the aforemen- 
tioned general structural phases [321 131] - The transition lines between the phases 
then depend indeed on both thickness and temperature. For low temperatures, 
the helical phase corresponds to polymers with low thickness, the sheet phase 
to a little higher thickness and the ring phase to the very thick polymers. 

Figure 8: Two N = 13 conformations with p = 0.73 (left) and p = 0.74 (right), 
which are not the ground-state conformations but have a just slightly higher 
energy than these. The conformations correspond to distinguished ground- 
state structures at shorter lengths (helical and crystallized on the sc lattice, 
cp. Fig. For reasons of better visibility, the thickness is not shown on scale. 



3.3 Selected Problems: Deeper Analysis and Remarks 
3.3.1 Simple cubic symmetry 

One of the basic observations within our polymer model relates to the fact, 
that the theory's parameter space of thickness p and r™ ln values possibly con- 
tains a likewise finite set of isolated and special parametric points, for which 

6 With "symmetric" we mean a somehow defined symmetry of torsion angles, e.g., torsion 
angles are pairwise equal relative to the middle of the chain. It is generally not essential for 
defining or distinguishing different "phases", but it is an interesting property and helps the 
understanding. Corresponding symmetry observables are not shown. 
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ground-state conformations exhibit a strict crystalline structure. These poly- 
mer conformations are characterized by a point set of monomer positions, that 
is frozen into a regular three-dimensional (3d) structure and none positional 
degree of freedom is left over. Additionally, these ground states inherit a fi- 
nite and possibly large ground-state degeneracy as there exist many ways to 
arrange the polymers sequence of monomers onto the frozen point set, without 
lifting the theory's energy to larger values. Typically, we then expect a number 
of ground-state configurations that increases exponentially, ng cx exp(ciV), with 
the chain length for long chains. Furthermore, the presence of crystalline ground 
states in the tube model possibly is attached to a triplet of global radius values 
r gc = 0.5774, l/y/2 = 0.7071, and r gc = 1, which denote the radii of circular 
polymers, that have an end-to-end distance of unity with exactly N = 3, 4, and 
N — 6 monomers on a circle. With numerical means it is then easy to show, 
that these particular radii result into triangular lattice (N = 3), simple square 
lattice (N = 4), and honeycomb lattice (N = 6) ground-state polymer point 
sets in two dimensions (2d), that is to say dimensional reduced tube polymer 
model. The phase space of 3d polymers as such is much larger than in 2d and a 
search at r gc = 0.5774 and r gc = 1 does not reveal any crystals for 3d polymers 
that would persist in the thermodynamic limit. However, at the particular value 
of r gc = 1/V2 and for = 1, we find stable 3d ground-state conformations 
with point sets of simple cubic symmetry and with large ground-state degener- 
acy. These crystals are likely to extend in the thermodynamic limit for large N 
values. 

The numerical simulations in case of the 3d thick tube model have been 
performed on chain length values N = 8, 9, ... , 32, 36, 40, 44, 48, 52, 56, 60, and 
N = 64. The global radius parameter was chosen to be p = l/v2 with the value 
r mm _ j f or £i ic p OS ition of the Lennard- Jones potential minimum^ We em- 
ployed parallel tempering simulations in the temperature interval 0.01 < T < 0.5 
with a temperature partition, that ensures acceptance rates around 0.5 for par- 
allel tempering swaps in-between neighbors in the temperature. For complete 
temperature interval coverage on a N = 32 chain a total of 39 temperature 
replica was needed. A single Monte Carlo run consists out of 10 9 x N monomer 
positional updates and from the ensemble of configurations the minimum-energy 
configuration was stored. A sequence of about 10 continuation runs for each 
chain length N with identical run parameters but with continued start configu- 
rations then yields an ensemble of about 10 ground-state estimates, and also the 
global ground-state estimate of the simulation. Finally, an adapted conjugate- 
gradient method was applied for refinement. We found that the efficiency of 
the Monte Carlo simulation in an attempt to populate statistically independent 
ground states rapidly degrades for chain length values N > 36 and therefore 
the longest chains (except the one at N = 36) are excluded from the further 
analysis. 

7 We also modified a single interaction term of the Lennard— Jones interaction in-between 
the polymers end to end: V = oo for r < 1 and V = for r > 1. This facilitates a perfect 
arrangement of the polymers monomers at the start and the end on a simple cubic lattice, if 
they prefer to be direct neighbors in space. 
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For a crystalline polymer conformation with simple cubic symmetry one can 
find a set of transformations, i.e., translations and orthochroncous rotations, 
that map the polymers point set to some point set of a simple hypercubic lattice. 
Denoting by x ? sc = mfe a with i = 1, . . , , N, a = 1, 2, 3, and with mf integer, a 
point set on a simple cubic lattice, the mean squared distance to a simple cubic 
lattice 

1 N 

4 = ^E( f r lymor -^ c ) 2 (7) 

i=l 

can be transformed to d1 c = 0, if also an appropriate set of mf values is chosen. 
Our main numerical result consists in the finding that all of our ground-state 
polymers at TV = 8,9 and for TV = 11, 12, ...,36 fulfill the numerically de- 
termined inequality d^ c < 0.000053 and therefore we observe a blatant simple 
hypercubic symmetry in the ground state of the theory for the considered chain 
length values to a high degree of numerical precision. A particular impressive 
ground-state conformation is displayed in Fig. |H1 where the TV = 27 polymer 
folds onto a 3 x 3 x 3 cuboid: Cub(3 x 3 x 3). For the TV = 27 chain, we per- 
formed a total of eight different continuation runs, which all yielded simple cubic 
symmetry in the ground-state estimates with values dg C < 0.000053. From these 
eight configurations, six were found with point sets isomorphic to Cub(3 x 3 x 3), 
however with five different mappings of the polymer sequence to the cuboid and 
with almost degenerate energy close to the ground-state energy. Given the nu- 
merical ability of the algorithm to identify different ground-state and near-by 
ground-state conformations, it appears unlikely that the true ground state has 
not been identified for the TV = 27 chain. A similar remark applies to all shorter 
chains. For purposes of future reference, and as a yard stick of our numerical 
precision we display in Table [1] ground-state energy density values eo = E/N 
from numerical simulations (second row) as a function of the chain length TV, as 
well as exactly calculated energy densities for various cuboids. It is noteworthy 
that ground states at TV = 8, 12 and TV = 18 also exhibit cuboidal point sets: 
Cub(2 x 2 x 2), Cub(2 x 2 x 3) and Cub(2 x 3 x 3) respectively. For the TV = 36 
chain the exactly calculated energy density of the Cub(3 x 3 x 4) cuboid under- 
shoots the numerical value significantly and, in fact, the cuboid was not found in 




Figure 9: A ground-state conformation for the TV = 27 polymer with p = 1/V2 
and rfj m = 1 (i.e., a ^ 1). The left picture shows the conformation in the proper 
scale, the right side shows the same conformation without proper thickness. 
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Table 1: Ground-state energy density values eo = E/N from numerical sim- 
ulations {second row) as a function of the chain length N, as well as exactly 
calculated energy densities eg (cuboids) for various cuboids (fourth row). 



N eo Manifold eo (cuboids) 



8 -1.0032 Cub(2 x 2 x 2) -1.0038 

9 -0.9670 

10 -1.0398 

11 -1.1052 

12 -1.2497 Cub(2 x 2 x 3) -1.2514 

13 -1.2263 

14 -1.2773 

15 -1.3179 

16 -1.4152 

17 -1.4451 

18 -1.5267 Cub(2 x 3 x 3) -1.5293 

19 -1.5165 

20 -1.5496 

21 -1.5729 

22 -1.6346 

23 -1.6612 

24 -1.7159 

25 -1.7465 

26 -1.7944 

27 -1.8385 Cub(3 x 3 x 3) -1.8433 

28 -1.8230 

29 -1.8282 

30 -1.8383 

31 -1.8750 

32 -1.8875 

36 -1.9224 Cub(3 x 3 x 4) -2.0036 



the numerical simulations. This again indicates the failure of our numerical algo- 
rithms for chains with length N > 36. Finally, it is also of interest to classify the 
secondary structures of compactificd ground-state conformations for the case of 
simple cubic symmetries. In particular we may consider U-turns (planar), and 
simple cubic helices (3-dimensional) [39], which both are chain segments of four 
monomers with bending angles of 90-degrees in-between consecutive segments 
on the simple cubic lattice. Using a pattern recognition program along the ideas 
of Tenenbaum ct al. [40j on a set of five different ground-state conformations 
for the N = 27 polymer we obtain rather low probabilities Pu-tum ~ 0.27 and 
-Psc-hdix ~ 0.22 for the occurrence of U-turns and sc-helices, respectively. 
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3.3.2 The a-helix region 



For the N = 8 and 9 polymer, we found a thickness region, where the a-helix 
is the ground-state conformation (see Sect. 13.2]) . Remember that we used the 
Lennard- Jones potential with a = 1 there, which sets the interaction length 
scale. There is nothing special with it, except that the potential just vanishes 
at the bond length, a fact that plays just a "second order" role, as we are not 
counting energy contribution from consecutive monomers at all. 

Because of the special role of the a-hclix in nature (besides its geometrical 
elegance), we will here try to track the a-helix not only in the thickness but also 
in the cr-direction of the "phase space" , i.e. we vary p and a independently in the 
vicinity of the assumed "a-rcgion" . Our results are displayed, exemplarily for 
N = 8, in Fig.[l0j We see, that the a-hclix occurs as ground-state conformation 
in a small, bounded thickness interval (0.66 < p < 0.75) right "before" an 
abrupt conformational change (depicted by the solid line) to cubelike structures. 
The transition line increases approximately linearly in the interaction-length- 
thickness plane, a dependence, which seems to hold generally for structural 
transitions in the vicinity. Following a perpendicular path, i.e., with increasing 
interaction lengths and decreasing thickness, the helices untwist smoothly. The 
dashed- dotted line together with the solid line define the region, where the a- 
helix is the ground state of the system (cp. insets in Fig. [I}. Note that for 
p < 0.66 and p > 0.75, a-helices are no ground states at all. 




0.65 0.67 0.69 0.71 0.73 0.75 0.77 

P 

Figure 10: N = 8: The a-p plane and ground-state conformations near the 
a-hclix. The left and right coordinates are connected just via r™ in = 2 1 / 6 cr. Sec 
text for details. 
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A further interesting transition is marked by the dashed line in Fig. [TUJ 
This line indicates the transition between the so-called simple cubic "lattice he- 
lices (ii)" and "(i)" |35], i.e., cuboidlike structures with parallel and antiparallel 
tails (remember that for r™ n = 1 and p = l/y/2, we observe the "crystalliza- 
tion" exactly at the simple cubic lattice, as mentioned in Sect. 13.2.11 and in 
detail discussed in Sect. I3.3.T1) . 

For the sake of completeness, the dotted line indicates a conformational 
change to some less interesting intermediate structure "between" a- and lattice 
helices and the arrows on the y-axes mark the line a = 1 investigated in the 
first part of this study. 

4 Summary 

The aim of this work was to take the simplest coarse-grained model for off- 
lattice polymers with explicit volume exclusion and to show, to which degree 
polymer crystallization can be understood even with this simplest model. We 
introduced the polymer volume using the concept of the radius of curvature 
which is indeed, in the first instance, a mathematical concept. In fact, it has 
been proven that "[it] is connected to various physically appealing properties of 
a curve. In particular, [it] provides a concise characterization of the thickness of 
a curve, [. . . ] as have been investigated within the context of DNA" [T^ , it was 
further successfully used in more complex models for proteins [TCI US OH] and 
it was finally shown, that this concept is effectively equal to a volume exclusion 
using two-point functions for polymers in good solvents [18] . 

Using sophisticated simulation techniques, we have analyzed systematically 
and in detail ground-state structures for the described model with fixed inter- 
action length. We have shown, for example, that already in this simplest model 
basic secondary structures like helices and sheets form. This statement is, due 
to the simplicity of the model, valid for various classes of polymers. Of course, it 
should be stated as well, that mentioned structures are not very stable against 
variations of the thickness, but this was not expected either. 

We investigated furthermore in detail the "neighborhood" of the a-helix by 
varying both, thickness and interaction length. Affirming above statement, it 
turned out that the a-helix exists as ground state only in a small, bounded 
area in the a-p space, but is surrounded by other helical and helicallike con- 
formations. 

It was of course known for a long time, that helices and sheets form within 
coarse-grained models including a somehow defined volume exclusion, but to 
our knowledge mainly for dedicated or less simple and not that general models. 
In some interesting works, for example, the strength of directionalized inter- 
actions pj] , explicit hydrogen bonds (4TJ [42] , solvent particles [20l [2TJ [22] , or 
the interplay between attractive interactions and packing |14| play a role. In 
particular the findings of [T3] also confirm the existence of not too long heli- 
cal structures in a specific homopolymcr model that is characterized by strong 
repulsive interactions between spheres. It is a common ansatz to investigate 
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and understand protein folding, stressing that we do not speak only of proteins 
but of a general class of polymers including proteins, at different abstraction 
(coarse grained) levels. It seems that at least parts of the general secondary 
structure formation can be attributed already to the simplest generic model 
for thick polymers. Obviously, these secondary structure segments have to be 
strengthened by further interactions in order to reach, for example, biologically 
relevant structure sizes, as tertiary effects set in at comparatively short chain 
lengths in simple coarse-grained models [T2J [2TJ [32] . 

The analysis of ground states is, of course, just a first step to an understand- 
ing of the model. In a subsequent work [32] we will, based on the knowledge 
of the ground-state conformations, focus on the thermodynamic behavior and 
conformational phases at finite temperatures. 
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